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Resumen 


Herrera-Díaz, I. E. Rodríguez-Cuevas, C., Couder- 
Castañeda, C., dx Gasca-Tirado, J. R. (enero-febrero, 2015). 
Modelación numérica hidrodinámico-hidrológica en zonas 
de inundación con presencia de infraestructura. Tecnología y 
Ciencias del Agua, 6(1), 139-152. 


Se presenta la modelación numérica computacional de 
la hidrodinámica superficial e hidrológica de la zona de 
estudio donde se pretende construir infraestructura para 
la exploración de hidrocarburos en las márgenes del río 
Grijalva, México; para ello se muestra el análisis de la 
información meteorológica, determinando los valores de 
intensidad de precipitación, temperaturas, evaporación 
y una estimación de los gastos, correspondiente al año 
2014; las zonas de inundación se estimaron con el uso de 
programas desarrollados en Fortran y Matlab, que resuelven 
las ecuaciones de Navier-Stokes-Reynolds para flujos a 
superficie libre y la intensidad de precipitación con la 
distribución de Gumbel, con parámetros estimados mediante 
el método de Momentos Ponderados con Probabilidad 
(MPP), con los cuales, mediante una malla numérica de la 
topografía, en conjunto con los valores meteorológicos, como 
condiciones iniciales y forzantes, se determina la magnitud 
de la inundación del área de estudio, así como la obtención 
de gastos, velocidades y el funcionamiento hidráulico de las 
obras de mitigación propuestas para preservar el balance 
hidrológico del sistema. 


Palabras clave: zonas de inundación, modelos hidrodinámi- 
cos, modelos hidrológicos, río Grijalva. 


Abstract 


Herrera-Díaz, 1. E., Rodríguez-Cuevas, C., Couder-Castañeda, 
C., 6 Gasca-Tirado, J. R. (January-February, 2015). Numerical 
Hydrodynamic-Hydrological Modeling in Flood Zones Containing 
Infrastructure. Water Technology and Sciences (in Spanish), 
6(1), 139-152. 


The computational numerical modeling ofthe surface hydrodynamics 
and hydrology of the study area is presented, where the construction 
of infrastructure to explore hydrocarbons is planned on the banks of 
the Grijalva River, Mexico. Meteorological information is analyzed 
to determine precipitation intensity, temperatures, evaporation and 
flows for the year 2014. Flood zones are estimated using programs 
developed in Fortran and Matlab to solve the Navier-Stokes-Reynolds 
equation for free surface flows. The intensity of precipitation was 
obtained with the Gumbel distribution with parameters estimated 
using probability weighted moments (PWM). Based on this, the 
magnitude of the flooding in the study area was obtained according 
to a numerical grid of the topography along with meteorological 
values such as initial and forcing conditions. Flows, velocities and 
the hydraulic functioning of mitigation works proposed to preserve 
the hydrological balance of the system were also obtained. 


Keywords: Flood zones, hydrodynamic models, hydrological 
models, Grijalva River. 
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Introducción 


La exploración de hidrocarburos en México es 
una de las principales labores de la industria 
de los energéticos. Contar con estudios que 
proporcionen información veraz y oportuna 
es primordial para una acertada toma de 
decisiones, 
actualizado sobre aspectos técnicos y socioe- 


ya que genera un panorama 


conómicos que influyen y coinciden en un 
mismo espacio-tiempo. 

Una de las herramientas que están contri- 
buyendo a los estudios previos hidrológicos 
son los modelos numéricos computacionales, 
aplicando la dinámica de fluidos (CFD, por 
sus siglas en inglés), el cual constituye una 
aproximación en el estudio de las ecuaciones 
de la dinámica de fluidos que se utilizan 
para describir diferentes fenómenos físicos 
relacionados con el movimiento de los fluidos, 
como los producidos por los flujos compresibles 
e incompresibles. 


Este trabajo radica en acoplar un modelo 


numérico hidrodinámico y un modelo 
numérico hidrológico para determinar el 
funcionamiento hidrodinámico de la zona 
propuesta para la colocación de pozos de 
exploración o explotación, así como de la 
infraestructura de proyecto (caminos, accesos, 
puentes, etc.); se tiene conocimiento previo de 
que la zona se comporta como una llanura de 
inundación por sus características topográficas, 
hidráulicas y geohidrológicas, por lo que se 
plantea como objetivo secundario analizar 
medidas de mitigación, para no alterar de 
forma sustancial el balance hidrodinámico de 


la zona seleccionada. 


Localización y características de la zona 
de estudio 


El sitio de estudio se encuentra en la margen 
izquierda del río Chilapa o Grijalva en el estado 
de Tabasco, México (figura 1), donde está 
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Coordenadas: 


a) X = 534 741.66, Y =2 025 068.62 
b) X = 534 840.18, Y =2 025 069.98 
c) X =536 188.14, Y =2 023 546.00 
d) X =536 549.00 Y =2 021 843.93 
e) X =535 051.37, Y =2 020 875.74 
f) X = 534 287.24, Y =2 021 392.14 


Figura 1. Localización de la zona de estudio (municipio de Centla, estado de Tabasco, México). 
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limitada por un polígono formado por canales 
en la parte inferior y superior, así como del río 
Tabasquillo en el extremo izquierdo. Existe una 
Obra civil, la cual se ubica en la parte central de 
la planicie que forma el sistema y es el objetivo 
principal del estudio. 

El trazo del camino propuesto hacia las 
instalaciones proyectadas del pozo se presenta 
en la figura 2. En el camino, al ser nuevo en 
el sistema ambiental, se contempla el uso de 
pasos de agua o alcantarillas para permitir el 
flujo hidrológico del sistema sin que presente 
alteraciones sustanciales a su estado original. 
La ubicación y el dimensionamiento de éstos 
se determinan en función de las magnitudes 
de velocidades calculadas y de los gastos 
provenientes de los resultados de la simulación 
y deben garantizar que el flujo hidrológico 
no se interrumpa y siga su cauce natural, 
a fin de dar cumplimiento a la NOM-022- 
SEMARNAT-2013. 

Desde el punto de vista hidráulico, se 
considera que la morfología general de la zona 
es sensiblemente constante en el tiempo; esto 
es, la zona de inundación permanece de forma 
continua cubierta porvegetaciónensutotalidad, 
por lo que las magnitudes de velocidades del 
flujo superficial no son considerables como 
para alterar la morfología del terreno y, por 
consecuencia, ésta se considera invariable para 
efectos de simulación. 

La información para la precipitación se 
recopiló de las cartas de isoyetas (SCT-DGST, 
2013) correspondientes a periodos de retorno, 
intensidad de lluvia y duración de la misma, 
para poder modelar el flujo superficial en 
diferentes periodos de retorno y de intensidad 
de lluvia; la información se condensó y se 
aplicaron aquellas lluvias con duración de 30 
minutos para distintos periodos de retorno 
(figuras 3a a 3e). 

Se considera un coeficiente de escurrimiento 
que varía de 15 a 25%, pues presenta una 
permeabilidad media y baja, en función del 
material que está constituido el suelo, además de 
considerar una vegetación de densidad media 
en promedio para todo el sistema ambiental. 


Coordenadas pozo 
Y =2 022 547 840 
X =535 017 000 
Geográficas 
18% 17' 35.951” N 
92% 40' 7.235” W 


UTM 


Río Grijalva 
km 3 + 640 


Alcantarilla tipo 1 km 1 + 300 


km 0+000 A 


Planta Esc 1: 7 500 


Figura 2. Trazo propuesto de la infraestructura dentro del 
sistema ambiental. 


Se proponen dos escenarios a estudiar en 
un año: el primero para la época de secas y 
el segundo para la de lluvias (Yin, Yu, Zhane, 
Wang, é Xu, 2013). Los datos para cada 
escenario se presentan en los cuadros 1 y 2. 

Con la información obtenida se estiman los 
gastos del río Grijalva para el año 2014, del cual 
se proponen usar los gastos para la simulación 
hidrodinámica en 2D del tránsito de avenidas 
en las épocas de secas y lluvias, acoplado con 
un modelo hidrológico (Mani, Chatterjee, « 
Kumar, 2014; Zhang, Long, Xie, Zhu, £ Wang, 
2007). En la figura 4 se presentan las descargas 
del río. 
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Figura 3a. Isoyetas de intensidad de lluvia para un periodo de retorno de 10 años y duración de 30 minutos. 
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Figura 3b. Isoyetas de intensidad de lluvia para un periodo de retorno de 20 años y duración de 30 minutos. 
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Figura 3c. Isoyetas de intensidad de lluvia para un periodo de retorno de 25 años y duración de 30 minutos. 
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Figura 3d. Isoyetas de intensidad de lluvia para un periodo de retorno de 50 años y duración de 30 minutos. 
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Figura 3e. Isoyetas de intensidad de lluvia para un periodo de retorno de 100 años y duración de 30 minutos. 


Cuadro 1. Resumen de valores obtenidos de los datos para alimentar el modelo hidrodinámico para el año 2014. 


- E Precipitación media Temperatura media Evaporación media 
Epoca del año 
mensual mensual mensual 
Secas 156.79 mm 26.4 “C 148.66 mm 
Lluvias 206.15 mm 26.2 C 117.33 mm 


Cuadro 2. Valores de escurrimiento y permeabilidad en el sistema ambiental. 


Época del año 


Densidad de vegetación 


Permeabilidad del suelo 


Coeficiente de 
escurrimiento en la zona 


(%) 
Secas Media-baja Media De 15a25 
Lluvias Media Media Mayor que 25 
Metodología 


El modelo hidrodinámico empleado resuelve 
las ecuaciones de Navier-Stokes simplificadas 
para flujos de superficie libre en coordenadas 
cartesianas; con la hipótesis de presión 
hidrostática y considerando los postulados de 
Reynolds son (Broomans, 2003): 


0U dv 0U dv 
— VU+W—= 


0 [f",0az] (1) 
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Figura 4. Estimación del gasto del río Grijalva para 365 días del año 2014. 


donde v, es el coeficiente de viscosidad efectiva, 
obteniéndose mediante la suma del coeficiente 
de viscosidad turbulenta y molecular v,=0,+0, . 


Stansby (2003) propone el siguiente modelo 
para resolver la viscosidad turbulenta: 


2 2 2 
D,= e 9 e +2 dv pr du , du 
Ox dy ox Oy 
1/2 


dar. 
OZ OZ 


donde la longitud vertical se define por 
(, = K(z — z,) para (z — z,)/8<2/K y £, = 
para A/K < (2 — z,)/8 < 1. K es la constante de 
von Karman con valor típico de 0.41; (z — z,) es 


la distancia desde el fondo; d es el espesor de 
la capa límite, y A es una constante con valor 
de 0.09. Las condiciones de frontera para la 
superficie libre y de fondo se representan por: 


fondo du $ V 7 + p? 
: e A (u), 
fondo 
ondo 00 e +o 
a 


donde Cz es el coeficiente de fricción de Chezy. 
Las componentes de la velocidad toman los 
valores de la capa adyacente a la interface 
agua-sedimento, definidas como: 


roo) P.. 
sup _ _ VFaire 
E te dz Asiento o Xx O, , 
z superficie Pagua 
sup dv Paire ] [ 
y e viento y y 
0Z superficie Pagua (4) 


donde p,.. = 1.29 kg/m” w, y 0, son 
las componentes horizontales en x y y, 
respectivamente, de la velocidad del viento a 10 
m de altitud. El coeficiente a,... (adimensional) 
se puede obtener utilizando la ecuación dada 
por Flather (1976) (UNESCO, 1981): 


=0.565x10"si]g|<5% 
Ss 


a. 
viento 


doo = (-0.12+0.137|0 |) 10? si 5 s [[)| 19.22 m/s 


viento 


a 


viento 


=2.513x 107 si p | >19.22 m/s (5) 
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La amenaza por inundación está en 
función de la ocurrencia del fenómeno y de su 
intensidad (Gurusamy éz Jayaraman, 2012). La 
intensidad, a su vez, se define en función de la 
velocidad y profundidad del agua, así como 
de la duración de la inundación. Por lo tanto, 
la definición de amenaza de inundación, en 
función del daño potencial (Zhou et al., 2014), 
toma en cuenta la probabilidad o frecuencia 
de ocurrencia de la inundación, así como los 
niveles o altura de la superficie libre del agua. 

El modelo hidrológico propuesto tiene como 
principio la relación precipitación-escorrentía e 
incorpora la generación de lluvias convectivas 
estocásticas (Gubarevat, 2011), que permite 
definir escenarios de amenazas específicos 
para la posterior evaluación del peligro de 
inundación asociado. Se emplea la distribución 
de Gumbel con parámetros estimados con 
el método de Momentos Ponderados con 
Probabilidad (MPP), considerado la mejor 
distribución que se ajusta a precipitaciones 
extremas para este caso (Jean, Kim, € Kim, 
2011): 


214) 


1(t,T) =at' (7) 


donde P, es la precipitación media; A es el 
área limitada por las poligonales de Thiessen o 
curvas de isoyetas; P,es la precipitación medida 


Ms 


en estación climatológica “1” o precipitación 
promedio entre dos curvas isoyetas; I es la 
intensidad de lluvia en mm/h. El inverso 
del periodo de retorno es la probabilidad de 
ocurrencia, es decir, la probabilidad de que se 
presente la intensidad de tormenta en un año. 
Los parámetros a, b, m, n son función de 
las características meteorológica de la zona. El 
hietograma acumulado está descrito como: 


H(t)=1(t,T)*t=af" (8) 


y la intensidad instantánea asociada como: 


dH(t) 
I(+) == (b+1)=at' (9) 
0-20.) 
Para completar la información faltante 
de registros pluviométricos, se emplea el del 
National Weather Service (Patro, Chaterjee, 


Mohanty, Singh, é Raghuwanshi, 2009), 
basado en la siguiente correlación: 
» 2217] 
N [1/7] (10) 


donde P es la precipitación que se extrapola, P, 
las precipitaciones de las estaciones existentes 
en el entorno de la que desconocemos, y r, es las 
distancias desde esta estación a las estaciones 
próximas. 

Greenwood, Landwehr, Matalas y Wallis 
(1979) introdujeron el método MPP como 
un procedimiento de estimación robusto, 
alternativo a los métodos de momentos y 
máxima verosimilitud. La estimación de 
parámetros por el MPP (Caeiro, Gomes, éz 
Vandewalle, 2014) comienza con los momentos 
M 


LjK 


M EXF (1-7) (11) 


LK" 

Siendo j y kenteros no negativos, el momento 
ponderado con probabilidad de orden (/, j, k) es 
proporcional al ¡-ésimo momento con respecto 
al origen de la estadística de orden (¡ + 1) para 


una muestra de tamaño n=k+j+1. 
La función de distribución de Gumbel está 


e 


donde a y m son los parámetros de la 


dada por: 


a a 


flop En 


Jess 


distribución, y x es la magnitud del evento 
extremo. 

De acuerdo con Greenwood et al. (1979), 
la expresión para la magnitud de la variable 
con un periodo de retorno T, según esta 
distribución es: 
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X, = m-aln |-In (1-+) (13) 
T 
donde: 
m=M, - ea (14) 
a Mo-2M, (15) 
In(2) 


Para e =0.5772 y M,=M, , , para k entero no 
negativo se estima como: 


1, LE) (16) 


Siendo n el número de datos históricos y 
x, el valor del dato 1, ordenados de mayor a 
menor. Operando la expresión para M, se tiene: 


M,=2Yx, (17) 
13 

y (18) 

4 tn-1 


con base en lo cual se determinan los valores de 
a y m con las ecuaciones (14) y (15). 


Generación de mallas numéricas 


Se elaboraron las mallas rectangulares con 
batimetría interpolada de la zona de estudio, 


con espaciamiento en ambas direcciones de 20 
metros (figura 5). 


Resultados de la simulación 


Detallando las etapas de simulación, éstas 
consisten principalmente en la modelación 
de dos escenarios en el año: la época de secas 
y la de lluvias; se escogen estas temporadas 
considerando que son las más representativas 
en el año. Las condiciones iniciales para cada 
temporada se expresan en el cuadro 3. 


Simulación del sistema ambiental en 
condiciones actuales (antes del proyecto) 
Los resultados de la simulación de la 
hidrodinámica superficial de los escenarios de 
secas con alturas de precipitación fue de 45, 80, 
120 y 156 mm, y de lluvias con alturas de 160, 
185, 197, 206 mm en la totalidad del sistema 
ambiental; las duraciones de lluvia fueron hasta 
de 2.5 horas continuas según registros. En la 
figura 6 se muestran las zonas potencialmente 
inundables. 

Se puede observar que las zonas potencial- 
mente inundables para una altura de precipi- 
tación de 45 mm presentan un porcentaje de 
inundación de 8% del total de la zona de estu- 
dio, mientras que para una altura de 206 mm, 
se Obtiene un porcentaje de área inundada de 
86%. 
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Num. de celdas dirección Y con incrementos de 20 m 
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Altura (m) 


20” A 
Núm. de celdas dirección X con incrementos de 20 m 


Figura 5. Malla numérica de la zona de estudio (izquierda, estado natural; derecha, con camino proyectado). 
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Cuadro 3. Resumen de condiciones iniciales y forzantes para la modelación en las diferentes etapas y escenarios propuestos. 


Precipitación o Evaporación Escurrimiento 
Etapa Temporada a Temperatura (*C) O GRAN 
Condición actual Secas De 45 a 156 26.4 148.66 0.020 
y c/ proyecto Lluvias De 160 a 206 26.2 117.33 0.037 
Et parida Permeabilidad Densidad de Gasto Grijalva Velocidad viento 
ps a (%) vegetación (m/s) laguna (m/s) 
320 mín 
E ia-baj 1 
Secas 30.00 Media-baja 1.070 med NE 4 
Condición actual A 
y c/ proyecto 1 900 máx 
Lluvias 20.00 Media 3 050 máx EAS 
hist 


Simulación del sistema ambiental con obras 
de proyecto 


Para esta simulación se contempla la construc- 
ción de un camino, considerando obras como 
alcantarillas y puentes, y/o pasos de agua; es- 
tas obras fueron ya propuestas y ubicadas de 
manera estratégica. Los siguientes resultados 
muestran el comportamiento del camino con 
las alcantarillas y, en especial, de la zona de 
inundación, donde se observa que el funciona- 
miento hidrológico del sistema ambiental no se 
ha alterado con estas medidas de mitigación. 

De los resultados de las zonas de inun- 
dación para los escenarios de intensidad de 
precipitación descritos en la figura 7, se aprecia 
que el nivel de la superficie libre originada por 
lluvia y por el desbordamiento del río Tabas- 
quillo no afectó el camino ni lo tapó. 

Para el escenario antes de proyecto, se 0b- 
servan las mismas zonas inundadas, así como 
algunas zonas con acumulación o represamien- 
to de agua. El porcentaje de zona inundable con 
el camino es de un 89%; es decir, 3% más que en 
condiciones normales antes de proyecto. 


Análisis hidrodinámico de las medidas de 
mitigación (alcantarillas y pasos) 


El camino, al ser nuevo, requiere de drenes 
o alcantarillas, para permitir el libre paso 
del flujo hidrológico en la zona de estudio. 
La propuesta establecida por el proyecto en 


cuanto al número, dimensiones aproximadas y 
ubicación de estas obras a lo largo del camino 
(figura 2) se resume a continuación: 

Se simuló el comportamiento hidráulico 
para los dos tipos de las alcantarillas descritas 
en el cuadro 4 bajo las condiciones de alturas 
de precipitación de 45, 80, 120, 156, 160, 185, 
197 y 206 mm, en la figura 8 se muestran los 
resultados para la simulación de 206 mm 
para la alcantarilla tipo 2 (1.05 x 1.05 m), 
considerando ser la situación más desfavorable 
y la que presenta mayor interés por las posibles 
magnitudes de velocidad que alcanza en la 
sección de la alcantarilla. La simulación del 
flujo a través de la obra de paso tiene como 
objetivo revisar los siguientes criterios: 


+ El área hidráulica debe garantizar el 


funcionamiento hidráulico las 
intensidades de lluvia más críticas, según 


registros y tiempos de retorno del sistema 


para 


ambiental. 

+ La densidad de la vegetación del sitio para 
cada alcantarilla será diseñada dando un 
cierto porcentaje por crecimiento de la 
vegetación (obstrucción del área hidráulica 
de un 65 a un 70%), por lo que se propone 
de un 15 a 20% más del área de diseño. 


El comportamiento hidráulico de la sección 
de la alcantarilla reporta velocidad promedio 
de 0.05 m/s, suficiente para permitir el paso 
del flujo, considerando un área hidráulica 
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a) Zona de inundación con altura de precipitación de 45 mm 
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g) Zona de inundación con altura de precipitación de 197 mm 
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b) Zona de inundación con altura de precipitación de 80 mm 
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f) Zona de inundación con altura de precipitación de 185 mm 


Altura (m) 
So 


180. Too 


60 


60 o 40 
O 20 
0.0 


1 , Núm. de celdas dirección Y con incrementos de 20m Núm. de celdas dirección X con incrementos de 20 m 
Núm. de celdas dirección X con incrementos de 20 m 


h) Zona de inundación con altura de precipitación de 206 mm 


Altura (m) 
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Figura 6. Zonas potencialmente inundables antes de proyecto. 


efectiva de 35%; esto es, el área restante que 
es obstruida se atribuye a la influencia de 
la vegetación que llega a cubrir la sección 
de dicha alcantarilla; esta condición se 
implementa en la simulación, considerando 
obstáculos en la sección que asemejen la 
condición de la vegetación. 

El valor de la amplitud de inundación en- 


tre la precipitación mínima y la máxima es de 


1.60 metros, el valor con el que se modeló la 
diferencia de niveles de superficie en los dos 
tipos de alcantarillas (de 0.91 y 1.05 m). Con es- 
tos resultados se estima que la magnitud de la 
velocidad es suficiente para generar circulación 
en la zona de inundación y brindar en menor 
cantidad alteraciones a la zona de estudio, pro- 
piciando condiciones similares a las originales 
antes del proyecto. 
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a) Zona de inundación con altura de precipitación de 45 mm b) Zona de inundación con altura de precipitación de 80 mm 
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c) Zona de inundación con altura de precipitación de 120 mm d) Zona de inundación con altura de precipitación de 156 mm 
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g) Zona de inundación con altura de precipitación de 197 mm bh). Zona deinundación conialtura de precipitación de 206 mm 
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Figura 7. Zonas potencialmente inundables, incluyendo el camino. 


De) 
a Cuadro 4. Ubicación y dimensiones de las alcantarillas en el camino proyectado. 

vu 

2 

o 

5 Di i Icantarill 

2 Kilometraje de ubicación Número de alcantarillas EAS cs 

jo) 

5 0+000 al 1+300 5 0.91 x 0.91 (tipo 1) 

EE 2+760 al 3+640 4 1.05 x 1.05 (tipo 2) 

E 

> y eo 

E Conclusiones esto se debe a que en la margen izquierda 
y del río, que es donde se encuentra la zona 
= 

= Los resultados muestran que el río Grijalva no de estudio, se localiza un bordo que protege 
E presenta una influencia superficial evidente; de alguna forma la crecida del río, con una 
3 


« vy80/0u72], 


Herrera-Díaz et al., Modelación numérica hidrodinámico-hidrológica en zonas de inundación con presencia de infraestructura 


Superficie del agua en secc. de 1.05; iteraciones = (491 000) 
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Figura 8. Comportamiento hidráulico de una alcantarilla de 1.05 x 1.05 m, con altura de precipitación de 206 mm. 


altura aproximada de 3.7 metros, y por lo 
tanto minimiza el riesgo de inundación. Sin 
embargo, el río Grijalva tiene una influencia 
sobre los canales superior e inferior, que con 
ayuda del río Tabasquillo permiten inundar 
la zona de estudio. Las zonas potencialmente 
inundables se consideraron de 0.5 a 3.5 metros 
de elevación, que ocupan un aproximado de 
80% del sistema ambiental; en estas zonas se 
proyecta localizar la infraestructura del pozo. 
En cuanto a la posible ubicación del pozo, 
éste se localiza en una zona baja, propicia 
a la colección de escurrimientos, donde 
la velocidad del flujo que se reporta en la 
modelación que pasa por las alcantarillas del 
camino alcanza magnitudes de hasta 0.08 m/s 
bajo condiciones de lluvias (intensidad de 
precipitación promedio mensual de 120 mm/h, 
con periodo de retorno de 100 años y duración 
de 30 min), por lo que la obra del camino y la 
pera no alteran sustancialmente la intensidad 
del flujo; en lo que se refiere a la dirección, al ser 
esta zona un bajo, la dirección no se modifica 
debido a que todo se concentra en el centro. 


El aspecto hidrológico primordial que se 
debe mantener es el paso del flujo superficial; 
las alcantarillas permiten este funcionamiento y 
sólo se debe vigilar que la vegetación no llegue 
a cubrir más de un 65% del área hidráulica 
propuesta; este criterio se plantea observando 
las condiciones actuales donde se ubicarán las 
alcantarillas en el camino. 
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